Estimating the final spin of a binary black hole coalescence 
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We present a straightforward approach for estimating the final black hole spin of a binary black 
hole coalescence with arbitrary initial masses and spins. Making some simple assumptions, we 
estimate the final angular momentum to be the sum of the individual spins plus the orbital angular 
momentum of a test particle orbiting at the last stable orbit around a Kerr black hole with a spin 
parameter of the final black hole. The formula we obtain is able to reproduce with reasonable 
accuracy the results from available numerical simulations, but, more importantly, it can be used to 
investigate what configurations might give rise to interesting dynamics. In particular, we discuss 
scenarios which might give rise to a "flip" in the direction of the total angular momentum of the 
system. By studying the dependence of the final spin upon the mass ratio and initial spins we find 
that our simple approach suggests that it is not possible to spin-up a black hole to extremal values 
through merger scenarios irrespective of the mass ratio of the objects involved. 

PACS numbers: 



I. INTRODUCTION 

Our understanding of the dynamics of a binary black 
hole system has advanced at an impressive pace in the 
last several years as numerical simulations have provided 
the missing piece for describing the complete dynamics 
of the system, starting from the late inspiral, and con- 
tinuing through the merger and ringdown of the final 
black hole. While these simulations have been exploring 
the parameter space of possible configurations of initial 
masses and spins of the black holes, the process is nec- 
essarily slow due to the time consuming nature of the 
simulations. This has led to a number of recent works 
aimed at producing recipes that allow for predicting, to 
a certain extent, what can be expected as the final state 
of the merged black hole from a given initial configura- 
tion. This is important, not only in order to decide which 
situations might give rise to the most interesting dynami- 
cal behavior, but also for the possibility of exploiting this 
information within a simulation. Two examples of this 
are (1) in the extraction of waveforms, where choosing 
an adequate background can reduce the errors in deter- 
mining observable effects, and (2) in helping to provide 
an analytic description of inspiral-merger-ringdown to be 
used in the construction of templates. 

To date efforts to produce such estimates have pro- 
gressed on several fronts. The effective-one-body (EOB) 
approach P, which maps the two-body problem to 
the dynamics of a test-particle in a suitably defined 
spacetime, provided estimates for the final black hole 
spin [j, that turn out to agree with current nu- 
merical simulations to within ^ 10%. More recently, by 
combining the EOB approach with test-mass limit pre- 
dictions for the energy released during the merger and 
ringdown phases, Ref. Q has refined these predictions 
for the non-spinning case to obtain final spins that agree 
to within ~ 2%. On another front, several fitting formu- 
lae for the final black hole spin have recently appeared in 



the literature. These formulae are based on the results 
of available numerical simulations along with a judicious 
exploitation of the natural symmetries in the problem 
and/or expectations from the test-mass limit (see e.g. 

Refs. Ii,i,[i3,[ni,[il,[il)- 

In this work we follow a different route based on first 
principles, although implicitly using numerical results 
along with intuitive arguments that during the merger 
and ringdown phases, the mass and angular momentum 
of the system are roughly conserved. Beyond this as- 
sumption, however, and in contrast to the approaches 
mentioned above, we make no explicit use of results 
from post-Newtonian approximations or numerical simu- 
lations. Our straightforward approach is based on simple 
estimates for the quasi-adiabatic inspiral, coupled with 
standard results for the angular momentum of a particle 
in a Kerr black hole spacetime corresponding to the final 
black hole. We obtain a closed form expression for the 
final spin of the black hole for arbitrary mass ratios and 
spins. We find that predictions from this simple-minded 
expression give results that agree reasonably well with 
the numerical simulations. This is yet further evidence 
of the rather simple behavior underlying the dynamics of 
binary black hole spacetimes. 

In spirit, our work is similar to that of Ref. 0, who 
used a point particle approximation on a Kerr back- 
ground to estimate the mass and spin of the final black 
hole. Our approach, however, contains an important dif- 
ference with that of Ref. , which allows us to make rea- 
sonable predictions even for comparable mass systems. 

Our approach and assumptions are described in Sec. [Til 
and we present a simple expression that estimates the 
spin of the final black hole. In Sec. IIIII we illustrate our 
results for several interesting scenarios, including spin- 
ning, precessing black hole binary systems, and compare 
with numerical simulations. We conclude in Sec. IIVI with 
final comments. 
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II. THE APPROACH 

We consider two widely separated black holes that 
can be well approximated by two Kerr black holes with 
masses and angular momentum parameters (m^, Oi). For 
simplicity we will first restrict our discussion to scenarios 
where the orbit stays within a plane (which we will refer 
to as the equatorial plane). In this case the orbital an- 
gular momentum and individual spins are orthogonal to 
the equatorial plane. Our goal is to obtain the angular 
momentum of the final black hole in terms of the initial 
configuration of the system by a phenomenological ap- 
proach rather than by evolving the system numerically. 
This will aid in finding particularly interesting cases that 
can be explored with numerical simulations. 

Achieving this goal certainly requires some compro- 
mises as the system undergoes a long dynamical process 
that comprises several stages: inspiral, merger and ring- 
down. While an accurate description of the whole process 
requires following the system completely, a back of the 
envelope estimate can be obtained by exploiting the fact 
that: (i) during the inspiral phase, the system evolves 
quasi-adiabatically, and (ii) during the merger and ring- 
down phases the total mass and angular momentum of 
the system change by only a small amount. These obser- 
vations are based on intuitive arguments but are strongly 
backed by post-Newtonian and perturbative calculations 
for the inspiral and ringdown, and by numerical simula- 
tions for the merger. 

We obtain a simple expression, based on first princi- 
ples, that can be employed to predict what the angular 
momentum parameter of the final black hole will be. This 
expression is obtained naturally when, inspired by the 
observations above, the following assumptions are made: 

• To first order, the mass of the system is con- 
served. Thus, the final black hole will have to- 
tal mass M = mi -\- m2. During the whole 
process the total radiated energy remains small 
(-Mradiatod^ 10%Mi„itiai); thus this assumptiou is 
justified to the level of approximation we seek. 

• The magnitude of the individual spins of the black 
holes will remain constant. Since both spin-spin 
and spin-orbit couplings are small, and radiation 
falling into the black holes affects the spins by a 
small amount, this is a safe assumption. Therefore 
the contribution to the final total angular momen- 
tum due to the individual black holes spins will be 
determined by the initial spins. 

• The system radiates much of its angular momen- 
tum in the long inspiral stage until it reaches the 
innermost stable circular orbit (ISCO), at which 
point the dynamics quickly leads to the merger of 
the black holes. During the merger the radiation 
of energy and angular momentum with respect to 
the mass and angular momentum of the system is 
small. Thus, to estimate the contribution of the 



orbital angular momentum to the final angular mo- 
mentum of the black hole, we adopt the orbital an- 
gular momentum of a test-particle orbiting at the 
ISCO of a Kerr black hole with a spin parameter 
of the final black hole. While adopting this value 
makes strict sense in the extreme mass ratio case, 
we will sec that it leads to predictions that agree 
quite well with the results of numerical simulations. 
That the point-particle approximation is able to 
capture key aspects of the two-bod y dy namics has 
also been observed in Refs. [1, [fl UU. [l4| when com- 
paring with results in Ref . . in Refs. [H, 
when comparing the close-limit approximation to 
numerical simulations, and more recently in stud- 
ies of unstable circular orbits in black hole merg- 
ers [13. 

Bringing all these assumptions together, we may express 
the angular momentum parameter of the final black hole 

a/ as. 
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where ioib is the orbital angular momentum of a particle 
at the ISCO of a Kerr black hole (with spin parameter 
a/)- 

Note that our assumptions differ in several ways from 
those of Ref. 0: (l)we keep the mass of the system 
constant, while Ref. Q adds a contribution to the final 
mass from the energy at the ISCO; (2) we keep the con- 
tributions from the spins of both bodies, while Ref. Q 
neglects the spin of the smaller black hole; and (3) we use 
the orbital angular momentum of the ISCO for a Kerr 
black hole with the final spin, while Ref. uses the ini- 
tial spin of the larger black hole. For extreme mass ratio 
cases, both approaches would give essentially the same 
answer; however our approach can be applied to general 
mass ratios and accounts for the orbital and both indivi- 
dal spin angular momenta when obtaining the final black 
hole spin. 

We can re-express our formula for a / in a more conve- 
nient form as 
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where Xi = ai/rm (e [—1,1]) and u ~ (mim2)/A/^ (e 
[0, 1/4]), and without loss of generality it is assumed that 
mi>m2. This equation provides a way to obtain a/ 
given rui and Xi- Since the expression for Lorb at the 
ISCO is known in closed form for equatorial orbits we 
concentrate first on such cases. This will allow us to 
investigate the viability of Eq. by comparing it to 
results obtained with numerical simulations. 

Adopting the expression for equatorial orbits in 
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FIG. 1: The final spin Uf/M vs u = mim2/M^ for initial 
black holes with extreme spin parameters x = 1- Clearly, as 
the mass-ratio approaches f = 1/4 (equal-mass case) the final 
expected spin parameter decreases. 
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where the upper signs correspond to prograde orbits and 
the lower signs to retrograde orbits. This function is to 
be evaluated at r = risco with 

nsco = A/ {3 + Z2 T [(3 - Zi)i3 + Zi + 2Z2)f'^] , 
where 
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The use of the prograde or retrograde case depends on 
whether the final spin is aligned or anti-aligned with the 
initial orbital angular momentum. Indeed one particu- 
larly interesting application of the above expression is to 
understand the direction of the final spin as a function 
of initial spins and the mass ratio. 



III. REPRESENTATIVE CASES 

We now concentrate on several representative cases 
that explore different interesting scenarios, and make 
contact with available numerical results. 



FIG. 2: Final spin o//Af vs u for initial black holes for sev- 
eral spin parameters x = 0,0.2,0.4,0.6,0.8,0.85,0.9,0.95,1. 
For spins > 0.948 the maximum final spin is achieved as the 
extreme mass ratio case is approached while for initial spins 
^ 0.948 the equal-mass case would give rise to the maximum 
final spin. 



A. Equal spin case 

A simple case that can be compared with existing sim- 
ulations is for equal spins (i.e. xi — X2 — x)- In this 
case equation ^ reduces to 



M 



A//2 



+ (l-2i.)x. 



(5) 



This equation allows us to determine the value of a/ as 
a function of v and x a-nd answer specific questions. Fig- 
ure [1] illustrates the behavior of the final spin parameter 
as a function of mass ratio when the individual spins of 
the initial black holes are maximal. The largest spin for 
the final black hole is achieved for the extreme mass ra- 
tio case. This coincides with the intuitive picture that a 
particle falling into an extreme black hole will have neg- 
ligible effect on the final spin, while a head-on collision 
of equal-mass extreme Kerr black holes will give rise to 
a final black hole with af/M = 1/2. 

This behavior, however, varies when considering initial 
spins less than maximal. For initially non-spinning black 
holes, intuitively the final black hole will also be essen- 
tially non-spinning for the extreme mass ratio case while 
it would have a non-trivial final spin for the equal-mass 
case. Figure [5] illustrates the spin of the final black hole 
as a function of the mass ratio for different values of x for 
spins that are aligned with the orbital angular momen- 
tum. We see that there is a critical value for the initial 
spins, X = 0.948, above which, the maximum final spin 
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FIG. 3: Final spin af/M vs v for initial black holes with 
initially non-spinning black holes (x = 0) as predicted by our 
simple model (filled squares), and as obtained numerically 
in Ref. Q (red diamonds) , by the EOB model combined with 
non-spinning test-particle predictions in Ref. 0] (green circle) 
and by numerical relativity combined with non-spinning test- 
particle expectations in Ref. [TO] (blue triangle). The largest 
final black hole spin is obtained in the equal-mass case. 



will increase as the mass ratio q ~ mxjmi > 1 increases. 
Below the critical value, the final spin will increase as 
the equal-mass limit is approached. Finally, at the criti- 
cal value any merger will leave the final spin essentially 
unchanged irrespective of the mass ratio. 

The case of black holes which are initially non-spinning 
can be compared directly with a number of simula- 
tions [m, HI m, m. Hi]. For cqual masses, the value 
predicted by Eq. Q is af/M = 0.66, which is quite close 
to the value af/M — 0.68 obtained by simulations of 
equal-mass, non-spinning black holes. Figurc[3]illustrates 
our predicted values for the final spin as the mass ratio 
is varied from the extreme mass ratio to the equal-mass 
case, along with the results from Ref. Excellent 
agreement is found with results from numerical simula- 
tions. Notice that the final black hole with the highest 
spin occurs in the equal-mass case as expected since the 
orbital angular momentum is maximized in that case. 

Another case that can be compared to simulations is 
for equal masses where the initial spins arc cither aligned 
or anti-aligncd with respect to the orbital angular mo- 
mentum. While simulations for close to maximally spin- 
ning black holes are a challenge, robust results exist 
for Xi G [—0.7,0.7]. Figure m shows the final value for 
the spin parameter predicted by Eq. ([5]) as x is var- 
ied, along with the values available from existing simula- 
tions. Clearly, a quite reasonable agreement is found for 
Xi € [—0.7,0.7]. Furthermore, we obtain for the extreme 




FIG. 4: Final spin af/M vs x for equal-mass black holes 
as predicted by our expression (filled squares), numerically 
obtained (filled diamonds) in Ref. [l2[ and from the EOB 
model (filled triangle) used in Ref. [3, but without taking 
into account the angular momentum released during merger- 
ringdown (thus, the latter curve overestimates the results and 
should be considered as an upper limit). Irrespective of the 
alignment or anti-alignment of the individual spins, the fi- 
nal black hole spin is aligned with the initial orbital angular 
momentum. 



cases a/(x = -1) = 0.2909M and a/(x = 1) = 0.9591M, 
which are, respectively, 14% and 0.01% away from the 
values reported by the fit formulae employed in Ref. [l^] . 
Figure 0] also shows results from the EOB model used 
in Ref. [5|, where the final black hole spin is computed 
at the end of the EOB plunge, disregarding spin-spin 
corrections and without taking into account the angular 
momentum released during mcrger-ringdown ^; thus, the 
values should be considered as an upper limit ^. 

A final interesting example that we report here is where 
the individual spins are anti-aligned with the initial or- 
bital angular momentum and the spin of the final black 
hole is zero. This occurs in our model when the orbital 
angular momentum remaining at the ISCO is exactly 



^ Note that the (conservative) EOB Hamiltonian used in Ref. 
differs from the one used in Ref. Q]. Whereas Ref. d adds 
spin effects to the (rcsummod) non-spinning EOB Hamiltonian, 
which is a deformation of the Schwarzschild spacetime, Ref. U| 
also includes spin effects in the resummation, obtaining an EOB 
Hamiltonian which is a deformation of a Kerr spacetime. 

^ For example in the non-spinning case (xi = in Fig. \^ by in- 
cluding test-mass limit predictions for the angular momentum 
released during the merger-ringdown phases, Ref. Q has reduced 
the difference from ~ 10% to ~ 2%. 
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FIG. 5: Relation between the mass ratio v and the individual 
spins of the black holes to achieve a final non-spinning black 
hole. Only mass ratios 0.183 (g^3.15) will allow for the 
final spin to be anti-aligned with the initial orbital angular 
momentum direction if the individual black holes spin fast 
enough. 
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FIG. 6: Final spin aj/M vs a for equal-mass black holes 
with XI = 0.584 and X2 £ [—0.584, 0.584] as predicted by our 
simple model (filled squares) and as obtained numerically in 
Ref. [ll[ (red diamonds). 



the contribution of the individual spins is determined by 



counteracted by the individual black hole spins. This sce- 
nario determines the border between cases in which the 
direction of the final angular momentum (or final spin) 
is determined by the initial orbital angular momentum, 
or by the direction of the initial black hole spins. In the 
latter cases the direction of the total angular momen- 
tum will "flip" , and inertial frames will see the direction 
they are dragged reverse as the system goes through the 
merger. This should be reflected in the waveforms, which 
would likely display an interesting behavior [2^ . FigurelS] 
shows the mass ratios and initial value of the individual 
anti-aligned spins required to yield a final non-spinning 
black hole. For a realizable situation of = CSm^ the 
mass ratio adopted should be 4 : 1, while for maximally 
spinning holes, we would predict a final spin of zero for 
the mass ratio q ~ 3.15. Note that this value is smaller 
than that predicted by Ref. of g ~ 4.23. By our 
arguments, any system with mass ratio q > 6.78 will 
undergo a fiip of the total angular momentum if the indi- 
vidual spins are equal and anti-aligned with the initial or- 
bital angular momentum and their spin parameters obey 
ai/rrii > 1/2. Notice however that the orbital separation 
at which this flip takes place also depends upon the mass 
ratio. This location can be estimated via simple Newto- 
nian arguments. A particle of reduced mass ^ = vM in 
a circular orbit about a central object of mass AI has an 
associated orbital angular momentum given by, 

^approx = y.VrM , (6) 
Thus, the distance at which Lapprox will be canceled by 



M2 



= (1-21^)1x1 



Thus, the following estimate 



r 

M 



(1 



(7) 



(8) 



can be used to determine the approximate distance at 
which the flip will occur. Notice in the limit — 
[q oo), that r — > oo so the flip has essentially 

taken place prior to any astrophysically interesting initial 
configuration. In such cases the final spin direction is 
determined by the spin of the large black hole. 



B. Unequal spin case 

We now illustrate the case with equal masses, but un- 
equal spins. Setting Xi = X ~ 0^X2 and ly = 1/4 in 
Eq. ^ yields 



Lorb X / -, I \ 

M M^+4(^ + ")- 



a/ 



(9) 



Figure [5] illustrates the value of the final spin parame- 
ter for equal-mass black holes and x = 0.584 and vary- 
ing a G [—1,1], while Fig. [7] illustrates the case with 
X = —0.584. In both cases we compare directly against 
the results of the simulations in Ref. [ll|. Our results 
differ at most by 8% with the reported results. We have 
also compared with simulations with mass ratio 3 : 2 
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FIG. 7: Final spin aj jM vs a for equal-mass black holes with 
Xi = -0.584 and X2 G [-0.584, 0.584] as predicted by our 
simple model (filled squares) and as obtained numerically in 
Ref. p4| (red diamonds). 



and either a\jm\ ~ —0.194 and aijmi — 0.201, or 
ai/mi = 0.193 and 02/7712 = —0.201 In these 

cases our simple model predicts a final spin of 0.617M 
and 0.671M, respectively, while Ref. [2^ obtains numer- 
ically a final spin of 0.640M and 0.704M. Finally, for 
equal-masses and ai/mi = —0.198 and 02/7712 = 0.198, 
Ref. [13, nil obtains a final spin of 0.697M, whereas we 
predict 0.663M. 



C. Generic spin configurations 

Until now our analysis has been restricted to cases 
where the orbital plane does not change in time. How- 
ever, we expect that the same arguments which lead us to 
Eq. ^ are applicable to more generic scenarios with pro- 
cessing orbits and arbitrary directions for the individual 
spin. A key difference in generic cases is that the or- 
bit at the ISCO will, in general, be inclined with respect 
to the final total angular momentum. For these cases 
the expression for the orbital contribution to the total 
angular momentum would require cither the numerical 
integration of generic geodesies in a Kerr spacetime or 
the use of the radial potential for quasi-adiabatic spheri- 
cal orbits [2^ . Alternatively, one can make use of the fit 
formulae presented in Ref. 0] to express Lor& in analytic 
form. 

The simplest possible extension of our method to more 
generic spin configurations can be formulated when the 
following assumptions (in addition to the assumptions we 



FIG. 8; At ISCO, the total angular momentum will be the 
sum of the orbital angular momentum and the spin vectors. 
The angle j3 is assumed to be given and remain fixed during 
the insprial prior to ISCO. The inclination angle t is solved 
for along with the magnitude of the total angular momentum. 



made in Sec. ^ are adopted: 

• The following quantities arc known 

{7771, 7772, 5*1, S'2,iorb} 

at some point of the inspiral (prior to the ISCO), 
where Lorb is a unit vector parallel to the orbital 
angular momentum. 

• Both the magnitude of the total spin Stot — Si + 
5*2, and the angle 6'ls between the total spin and 
the direction of the orbital angular momentum Lorb 
will remain constant up to the ISCO. 

We notice that in general 0ls and l^totl can change 
during the evolution so this assumption might not hold 
to a tolerable level. However, we know from the post- 
Newtonian spin-precession equations [sot that this as- 
sumption is valid in two cases, notably (i) equal-mass 
double-spin binary systems (when spin-spin terms are 
neglected) and (ii) unequal-mass single-spin binary sys- 
tems [H]. 

From the initial conditions, we compute 



'S'tot = l'5'totl, 

Lorh ■ >S'tot 



COS f LS 



\Su 



(10) 
(11) 
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In Fig. [H we show the total spin and orbital angular 
momentum at the ISCO, where Jf = Lorb + •S'tot- As 
before we can then obtain the final spin of the black hole 
as Hf = Jf/M. More explicitly, if we decompose the 
vectors along the directions parallel and orthogonal to 
the final spin, we obtain 

iorb(t, a/)cosi + S'tot cos(6'ls - ') = Maf, (12) 
iorb(i, a/) sint - S'tot sin (6'ls - ') = 0, (13) 

where Lorb = |£orb|- Eqs. (|12m3p can be solved to derive 
the magnitude of the final spin af and the inclination 
angle t at the ISCO. For simplicity, we can compute the 
orbital angular momentum of the inclined orbit using the 
fit formula of Ref . 0] ■ 

iorb(t,a/) = ^(1 + cosOLP'bKsco'a/) + 

i(l-cosO|i^?bWsco,«/)l, (14) 

where L^J^ and L™^ are given by Eq. ([3|) for prograde 
and retrograde orbits, respectively, and ''fsco ^'^'^ ''iSCO 
are the corresponding ISCOs. 

With the above procedure, we can compare the ob- 
tained estimates with available numerical results. In Ta- 
ble m we list the estimates from our approach together 
with the results obtained from numerical simulations in 
Rcfs. [1, [111, [13] where several equal-mass double-spin 
and a few unequal-mass single-spin precessing configura- 
tions have been studied. Note that in cases Ci and C3, 
the total spin is zero (i.e. the individual spins arc equal, 
but opposite). In any case where the total spin is zero, 
our approach will predict the same result (independent 
of the magnitude of the individual spins), which given as 
a function of mass ratio by the bottom curve in Fig. [21 
Note that the cases Ce, C7 and Cs all have the same Stot 
and thus our approach predicts the same af/M and 
L. Moreover, in cases Ci and C3, the masses are equal 
and the total spin is zero, but the individual spins are 
not on the orbital plane (corresponding to the large kick 
configurations in [2^, [s^l ) • In these cases our approach 
predicts the same a//M and l. Once again in these more 
complex physical scenarios the agreement of our simple 
approach is reasonably good. 

We find it important to stress again that it may be 
that for more general processional cases — in which both 
black holes carry spin and their masses are not equal, or 
when the system undergoes a transitional precession [26[ 
— , the above approach may not provide a good approx- 
imation. We plan to investigate more generic cases in 
more detail in the future. 



IV. CONCLUSIONS 

We have presented an approach to obtain a simple ex- 
pression to estimate the spin of the final black hole pro- 
duced through a merger of orbiting binary black holes. 



Case [Reference] 


ai/M 


t(deg) 


af/M 


t(deg) 


Ci [31] 


0.67 


- 


0.66 





C2 [3]J 


0.72 


- 


0.71 


23 


C3 [3]J 


0.68 


- 


0.66 





C4 [3]J 


0.73 


- 


0.71 


23 


Cs [31] 


0.64 


- 


0.61 


34 


Cs [31] 


0.81 


- 


0.82 


14 


C7 [3]J 


0.80 




0.82 


14 


Cg [3JJ 


0.80 




0.82 


14 


SP3 [8] 


0.72 


18 


0.70 


21 


SP4 [8] 


0.81 


10 


0.80 


13 


SPo [32] 


0.50 


33 


0.48 


35 



TABLE I: In the second and third columns we list results from 
numerical simulations of precessing, spinning black holes, for 
the magnitude of the final spin and the angle between the 
final spin and the direction of the orbital angular momentum 
(if available). We use data from Ref. [sH (we indicate with 
Ci, i = 1, ...8 the spin configurations from left to right in Table 
I of Ref. [3l|]) and from Refs. [1, 113 (see Table I in each case). 
In the last two columns we list the estimates obtained with 
our approach. 

In this work wc have concentrated on several especially 
interesting cases, but others can certainly be explored as 
well. 

Notice that our work is complementary to recent works 
aimed at giving fit formulae for different physical quan- 
tities based on the results of numerical simulations (see 
e.g., Refs. MM [m, [13, [il), and to the FOB predic- 
tions [3, [3j [1, [S|- Our expression, however, does not rely 
directly on the simulations, but rather on a simple ap- 
proach based on first principles. On the other hand, it 
has an inherent amount of error due to its simple assump- 
tions. Confronting our predictions with available results, 
we find that they agree rather well considering the limi- 
tations of our simplistic approach. This fact gives further 
evidence to the rather simple behavior describing the dy- 
namics of orbiting black holes. The expression presented 
in this work can be employed to predict the outcome of 
the simulation for a large number of cases not yet studied 
and can help determine which parameter choices might 
give the most interesting results. 

For example: 

• For individual spins aligned with the orbital an- 
gular momentum and Oi > 0.948mi the final black 
hole spin is larger as v decreases. However for indi- 
vidual spins Oi < 0.948TOj, as ly is increased the final 
spin will be greater. This transition number has an 
inherent error due to our approximations, and thus 
should not be considered sharp. Nevertheless, we 
would expect a transition to occur near this value. 
Hence our simple-minded model suggests that it is 
impossible to spin-up a black hole through mergers 
to its extremal value even in the ideal case where 
both spins and the orbital angular momentum are 
aligned, and no dissipative effects exist to reduce 
the final spin. As indicated by Fig. ^ the only 
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way to get a final maximal spin requires an already 
maximally spinning black hole merging in extreme 
mass-ratio situations. Any other alternative in the 
highly spinning cases (a^ > 0.9487ti,;) will cause the 
final spin to decrease. If (a^ < 0.948mi) the final 
spin will only increase up to a value a sa 0.95. Af- 
ter this state, any further merger will essentially 
leave this value unchanged^. Consequently binary 
black hole systems would not give rise to an orbital 
hang-up due to having J > after the ISCO. 

• The direction of the orbital angular momentum 
of a system with arbitrary spins (perpendicular 
to the orbital plane) determines the final spin for 
i/> 0.183. For 0.183 however the final spin di- 
rection will depend on how large the individual 
spins are; this can give rise to a final black hole 
whose spin opposes the initial orbital angular mo- 
mentum direction. This scenario should give rise 
to an interesting phenomenology in the resulting 
waveforms. We stress again that this critical value 
is approximate given our simple assumptions, but 
such a critical value must exist. 

As discussed in Sec. IIII CI our approach can be gen- 
eralized to spinning, processing binaries. As long as the 
total spin of the system and the angle between the total 
spin and the orbital angular momentum arc preserved 
during the evolution, we have been able to compare our 
results to numerical simulations of precessing binary sys- 
tems, obtaining good agreement. We plan to carry out a 
more thorough study of more generic spinning, precessing 
binaries in the future and investigate further extensions 



of our approach. 

Some other applications of this approach would be to 
employ the predicted values in order to adopt a rea- 
sonably close background for perturbative approaches 
to study the after-merger epoch or to compute physical 
quantities with respect to this background like gravita- 
tional radiation. Additionally it can be used to aid in 
providing the quasi-normal mode frequencies (which are 
a function of the final black hole mass and spin) to be 
used in analytically matching the inspiral to the merger 
to the ringdown. 
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